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Abstract 

We develop a scheme for building the scalar exchange-correlation (xc) kernel of time-dependent 
density functional theory (TDDFT) from the tensorial kernel of time-dependent current density 
functional theory (TDCDFT) and the Kohn-Sham current density response function. Resorting 
to the local approximation to the kernel of TDCDFT results in a nonlocal approximation to the 
kernel of TDDFT, which is free of the contradictions that plague the standard local density approx- 
imation (LDA) to TDDFT. As an application of this general scheme, we calculate the dynamical 
xc contribution to the stopping power of electron liquids for slow ions to find that our results 
are in considerably better agreement with experiment than those obtained using TDDFT in the 
conventional LDA. 

PACS numbers: 71.15.Mb 
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I. INTRODUCTION 

Starting with the pioneering work of Runge and Gross,- time-dependent density func- 
tional theory (TDDFT) has evolved into a powerful tool for studying excitations in atomic, 
molecular, and condensed-matter systems.-*^'^*^ In the linear-response regime, the key quan- 
tity of TDDFT is the dynamical exchange and correlation (xc) kernel /xc(r, r', defined as 
the Fourier transform with respect to time of the functional derivative 

/,,(r,r,t-t) = ^-^-^, 

where Vxc and n are the time-dependent xc potential and particle density, respectively. 
In contrast with the xc potential in static DFT,-i^ the dynamical xc potential is strongly 
nonlocal with respect to space coordinates,^ to the point that a local-density approximation 
(LDA), understood as the zeroth order term in a regular gradient expansion, does not exist. 
Indeed, the use of LDA to treat genuinely dynamical effects (i.e., effects not captured by 
the adiabatic approximation) is known to lead to severe contradictions within the theory.— 
Despite impressive successes of TDDFT, there is still the want of a scheme for including 
nonlocality in xc kernels, accurate enough and practically convenient in applications.- 

Contrary to ordinary TDDFT, the time-dependent current density functional theoryii 
(TDCDFT) is known to allow a consistent LDA, which is believed to be of about the same 
level of accuracy for the time-dependent phenomena as the standard LDA is for ground- 
state properties. In many concrete applications, however, (including the calculation of the 
stopping power of electron liquids described below) it is the scalar xc kernel of the ordinary 
TDDFT, rather than the tensorial xc kernel of the TDCDFT, that naturally enters the 
equations describing the many-body effects. 

In this paper we exploit the fact that TDDFT and TDCDFT would be completely equiv- 
alent if the exact xc functionals were known, to construct a nonlocal approximation for the 
scalar xc kernel of TDDFT starting from the LDA for the tensorial xc kernel of TDCDFT. As 
we shall show below, the resulting nonlocal xc kernel of TDDFT satisfies the exact zero-force 
sum-rule, the violation of which within LDA to TDDFT had once provided the motivation 
for introducing TDCDFT.— 

We believe that our new nonlocal xc kernel has a broad range of potential applications, 
particularly in transport theory. As a first demonstration of its usefulness, we present here 
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the results of calculations of the stopping power of an electron liquid for slow ions, wherein 
we find that the contribution of the many-body dynamical xc effects is not only numerically 
important, but also leads to better agreement with experiment when the new nonlocal 
expression for f^c is used in lieu of the conventional LDA. 

The organization of this paper is as follows. In Sec [ITl we derive a formula expressing the 
exact scalar xc kernel of TDDFT through the exact tensorial xc kernel of TDCDFT and the 
Kohn-Sham current density response function. In Sec. Illll we summarize the formal TDDFT 
of the stopping power of an electron liquid for slow ions and discuss the difficulties the LDA 
runs into. In Sec. HV] we give the details of our calculational procedure, present results and 
their discussion. Sec. |V] contains our conclusions. Appendix is devoted to the interrelations 
between scalar and tensorial zero-force sum rules within the exact and approximate theories. 

II. SCALAR XC KERNEL OF TDDFT FROM THE TENSORIAL XC KERNEL 
OF TDCDFT 

We start from the expression of the xc kernel of TDDFT-i^. (the dependence on r, r', and 
uj is implied) 

fxc = Xks - X"' - (1) 

where x is the longitudinal density response function and xks is its single-particle Kohn- 
Sham (KS) counterpart. Similarly,— 

fxc,ij = XKS,ij " Xij ~ ^^*|r _ r'l 

where fxc,ij is the tensorial xc kernel of the TDCDFT, Xij ^^id XKS,ij are the many-body 
current density response function and its single-particle KS counterpart, respectively. In- 
verting the relation between the tensorial current density and the scalar density response 
functions 

X = -;^V. ■ x^J ■ V;., (3) 

we can write 

X-^ = -^V-V•(LxL)"^VV-^ 
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where L is the longitudinal projector operator Lij = VjVjV ^. Using a simple operator 
identity 

where T = 1 — L is the transverse projector, we can write for the inverse scalar response 
function 

and similarly for Xks- Using Eqs. ([1]), ([2]), (jl]), and the KS counterpart of the latter, we 
readily arrive at 

X {xks-L) -Xks {txksT) • VV-^ (5) 

Equation ^ is our desired and central result: it expresses the scalar xc kernel of TDDFT 
in terms of its tensorial counterpart of TDCDFT. 

In the case of a hounded system, the exact scalar xc kernel satisfies the zero-force sum- 
rule— 

j /,,(r, r', u) V'no(r') dv' = VV,,{t), (6) 

where no(r) and K:c(r) are, respectively, the ground-state density and the xc potential. On 
the other hand, the tensorial f^c and xks satisfy the corresponding sum-rules of TDCDFT.— 
In Appendix |X] we prove an important result that with any approximation to the tensorial 
fxc satisfying the zero-force sum rule, the corresponding scalar xc kernel of Eq. ([5]) satisfies 
the sum-rule ()6]). 

III. TDDFT OF THE STOPPING POWER OF ELECTRON LIQUID 

We now illustrate the usefulness of Eq. (j5]) by applying it to the problem of the stopping 
power of an electron liquid for slow ions. 

A. Formal TDDFT of the stopping power of electron liquid for a slow ion 

The stopping power dE/dx is the loss of energy per unit path of an ion moving through 
the electron liquid. The constant of proportionahty between dE/dx and the ion velocity 
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(for low velocity) defines the friction coefficient Q, which can be written as— 

Q = Qi + Q2 

where Qi and Q2 are the single-particle and the many-body dynamical xc contribu- 
tions, respectively. The single-particle (binary-collisions) contribution Qi can be expressed 

Qi = nokFatrikp), (7) 

where kp is the Fermi wave-number, atrikp) is the transport cross-section of the elastic 
scattering in the KS potential of an electron at the Fermi level, and uq is the electron liquid 
density in the absence of the ion. As shown in Ref. [ij, keeping the Qi part of the friction 
coefficient only is equivalent to using the adiabatic version of TDDFT. 
The many-body dynamical xc contribution Q2 is given by^ 



Q2 = - j [Vrno(r) • v][Vr'no(r') ■ v 
9Im/^c(r,r',a;) 



dr dr', (8) 

a;=0 



duj 

where fxd'^i r', ^) is the scalar xc kernel of the inhomogeneous many-body system of an ion 
at rest in electron liquid and v is the unit vector in the direction of the ion velocity. 
In the following, we focus on the calculation of Q2- 

B. Contradiction inherent in the LDA 

The simplest approximation, namely the LDA to the ordinary TDDFT,— amounts to 
setting 

/.,(r,r',.;) = /,\^[no(r),^]5(r-r'), (9) 

where /^^^(n, tu) is the g — limit of the longitudinal xc kernel of a homogeneous electron 
liquid of density n. By spherical symmetry, substitution of Eq. ([9]) into Eq. ([8]) yields 

00 

Q2 = -— / dr[rno{r)] '- . (10) 

6 J OOJ u}=0 
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In the limit of zero density of the electron liquid uq ^ 0, the independent-electron part 
Qi of Eq. ([7]) vanishes, but Q2 of Eq. ffTOj) gives a finite value, because the gradient of 
the ground-state density no(r) of an isolated atom is not zero and 9Im/^^^(n, ci;)/9co'|^=o is 
negativeA^ Thus, LDA to the scalar /^^ yields a finite friction coefficient even in the absence 
of the electron gas, indicating an obvious flaw of the approximation. Table [I] shows this error 
quantitatively for a number of atoms in comparison with friction coefficient at = 2.2. 

Atom He Be C O Ne Mg Si 
Q(r, = 00) 0.04 0.11 0.17 0.24 0.30 0.36 0.43 
Q(r^=2.2) 0.34 0.43 0.70 0.46 0.16 0.15 0.54 
% 12 25 24 52 188 240 80 

TABLE I: Inaccuracy of LDA to TDDFT: Friction coefficient of free space (r^ = 00) and that of 
an electron Uquid of = 2.2 for several atoms. Line 3 is the ratio of lines 1 and 2 (%). 

To check that Eq. ([5]) resolves this problem of the finite friction coefficient of free space, it 
is sufficient to notice that an isolated atom is a bounded system and hence the sum-rule ([6]) 
holds, which, substituted into Eq. (IHl), yields zero identically. Moreover, LDA to TDCDFT 
satisfies the zero-force sum-rule by construction,— ensuring, as is shown in the Appendix, 
that Eq. ([6]) holds for a bounded system even if the local version of TDCDFT is used in 
Eq. ©.^ 



IV. CALCULATIONAL PROCEDURE, RESULTS, AND DISCUSSION 

Our numerical procedure is to evaluate and invert operators entering Eq. ([5]) on an ortho- 
normal set of radial basis functions. For XKS/iji we have employed the standard method of 
using the static KS orbitals to build the independent-electron response function 

XKS,ij{r,r\uj) = ^?2o(r)(5(r - r) 6ij - 

E [C(r)V.^,(r) - ^,(r)V.C(r)] 

^ uj ~ ep + e^ + tri 

X [rpir'W^i'.ir') - ^.(r')V;^;(r')] , (11) 
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where V'a(r) and ea are the single-particle wave-function and eigenenergy, respectively, in the 
state a, and fa is the occupation number of this state. For /xc,ij, we use LDA to TDCDFT 



as22 



ic 



where 



/xc,ifc(r,r',u;) jfc(r')cir' = — x 



(r,cu) = -eL[no(r)]Vu.(r) 



lUJ 



exc{n) is the xc energy density, 



(12) 



(13) 



(14) 



is the stress tensor, and u(r) = j(r)/'n,o(r) is the velocity field. The viscosity coefficients are 
given by^ 



n 



lix) 

n 



J X( 



lUJ 



(15) 



and fl^^rp[n,uj) is the transverse xc kernel of the homogeneous electron liquid of density n. 

In Fig. [H we plot the results at = 1.59 corresponding to the valence electron-density 
of carbon.30 It is instructive that within 1 < < 14 both TDCDFT and TDDFT give 
virtually the same result, which, we believe, is generally true for light atoms in high-density 
electron liquid. Then, at higher Zi, rather abruptly, the dynamical xc contribution almost 
vanishes in our present calculation, which can be understood qualitatively recalling that for 
heavy atoms in the electron liquid the charge-density distribution is close to that of isolated 
atoms, and hence Q2 should be small. 

In order to compare our result with experiment, we plot the friction coefficient at = 2.2 
versus the atomic number of a moving ion in Fig. [2l Results of the calculations with neglect 
of the dynamical xc [Eq. LDA to TDDFT [the sum of Eq. ^ and Eq. and 
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FIG. 1: (Color online) Friction coefficient of electron liquid of carbon density (r^ = 1.59) versus 
the atomic number of an ion. Green squares are results with neglect of the dynamical xc [Eq. ([7])]. 
Blue triangles are results with the dynamical xc included within LDA to the conventional TDDFT 
[Eq. (jlOp ]. Red circles are results with the dynamical xc included with use of Eqs. ^ and (|12p -(|15 p . 

LDA to TDCDFT [the sum of Eq. © and Eq. ([8]) with /^^ given by Eqs. ([5]) and (fT2|)- 



ITSll ] are shown, together with the experimental data of Ref. 121| for ions moving with the 
velocity of 0.5 a.u. at the distance of 1.2 a.u. from the last atomic plane of the (111) 
surface of aluminum. The inhomogeneity of the electron density an ion travels through is 
weak under these conditions, and we have used estimated experimentally.— Moreover, the 
experimental stopping power is predominantly electronic since the trajectory of an ion is well 
separated from the lattice atoms. All together, these two conditions justify the comparison 
with the theory within the electron liquid model. The non-monotonic dependence of the 
friction coefficient on the atomic number of the ion (so-called Zi-oscillations) is known to 
result within the single-particle theory from the competition between the increase in the 
electron liquid-ion interaction with the growing charge of the nucleus of the ion and its 
decrease due to the screening by shells of bounded electrons of the pseudo-atom as well as 
its resonant states.— 

While LDA to TDDFT (triangles in Fig. [2]) largely overestimates the friction coefficient 
at Zi > 5, the results using Eq. ([5]) (circles in Fig. [21) are in good agreement with the 
experiment in a wide range of 3 < < 18. A deviation occurs at small and large Zi, where 
the experimental friction coefficient is lower than the independent-electrons calculations 
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FIG. 2: (Color online) Friction coefficient of electron liquid of = 2.2 versus the atomic number of 
an ion. The green squares are the results with neglect of the dynamical xc as obtained from Eq. ([7]). 
Blue triangles are the results with the dynamical xc included within LDA to the conventional 
TDDFT as obtained from Eq. pup . Red circles are results with the dynamical xc included with 
use of Eq. ([5]) for xc kernel and Eqs. (IT2] )- ()15p of LDA to TDCDFT. Black solid squares are the 

211 for ions {v = 0.5 a.u.) moving at a distance of 1.2 a.u. 



measured stopping power of Al of Ref. 



from the last atomic plane of the Al (111) surface. 

(open squares in Fig. [2]). This feature has recently been reported as due to the finite velocity 
of ions.— Hence it is an effect of the deviation from linear dependence of stopping power 
on velocity. The same effect gives a positive contribution at 8 < < 12, suggesting that 
combined with the many-body effects of the present theory the agreement with experiment 
can be further improved.^ In the range 13 < Zi < 17, the dynamical many-body effects 
seem to be solely responsible for the enhancement of the friction coefficient compared with 
the independent-electron theory. 

In Fig. [31 we plot the friction coefficient at = 2 versus the atomic number of ions in 
the range 5 < Zi < 39. This is compared with the available measured stopping power for 
ions with the velocity of 0.68 a.u. channeled along the (110) direction in gold. Due to the 
channeling, collisions with the lattice atoms again do not give significant contribution to the 
stopping power. It must be noted, however, that under channeling conditions the assumption 
of the nearly constant electron density is an uncontrolled approximation. One important 
qualitative conclusion we can draw from comparison of the theory and experiment in this 
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FIG. 3: (Color online) Friction coefficient of electron liquid of = 2 versus the atomic number of 
an ion. Green squares are the results with neglect of the dynamical xc [Eq. d?])]. Blue triangles 
are the results with the dynamical xc included within LDA to the conventional TDDFT [Eq. (|1U|) ]. 
Red circles are results with the dynamical xc included with use of Eqs. ([5|) and ()12p - (jl5p . Black 



solid squares are the measurements from Ref. 



231 of stopping power of Au for ions {v = 0.68 a.u.) 



channeled along the (110) direction. The black dotted line is the calculation of Ref. [2^ with the 
dynamical xc included within the linear-response theory of the homogeneous electron gas. 

case is that for Zi > 22 the role of dynamical xc effects becomes negligible in both experiment 
and the present theory, while LDA to the TDDFT yields these effects largely overestimated. 
Similar to the dip in Fig. [2], the underestimated theoretical values at 7 < Zi < 12 can be 
attributed to the effect of finite velocity.-^ However, within the range IQ < Zi < 19 the 
dynamical xc contribution is too small to account for the onset at the experimental data, 
nor can the persistent enhancement of the friction coefficient in this range be attributed 
to the effect of finite velocity within the independent-particle theory. Further studies are 
required to elucidate the nature of this onset, the inhomogeneity of electron density being 
the most plausible cause. The dotted line in Fig. [3] represents the friction coefficient of silver 
obtained with the dynamical electron-electron interactions included in Ref. 
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within the 



framework of the linear-response theory of the homogeneous electron gas. In the case of the 



non -deg enerate plasma, an approach similar to that of Ref. l25| has been reported in Refs. 
and 27. 
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V. CONCLUSIONS 



We have rigorously expressed the dynamical xc kernel fxc{T^, r', uj) of TDDFT in the terms 
of its TDCDFT tensorial counterpart and the Kohn-Sham current density response function 
of independent electrons. Then, using the local density approximation to TDCDFT, we have 
built a nonlocal approximation to fxc{i^,r' ,uj) which satisfies the exact zero-force sum-rule 
for bounded systems. We believe that our new approximation will be broadly applicable to 
a variety of problems in electronic transport theory. 

As a first application, we have calculated the dynamical xc contribution to the stopping 
power of an electron liquid for slow ions. In doing so we have resolved a basic difficulty of 
the conventional LDA - the finite friction coefficient of free space - and we have improved 
the overall agreement between theory and experiment. 

APPENDIX A: RELATION BETWEEN SCALAR AND TENSORIAL ZERO- 
FORCE SUM-RULES 

The following sum-rules hold^ for the exact tensorial xc kernel 



and for the KS and the interacting current density response functions, respectively. 



where Vo(r) is the bare potential. 

In this Appendix, we prove that for any approximation to fxc satisfying the sum-rule 
( lAll) . the corresponding scalar f^c obtained through Eq. ([5]) satisfies the sum-rule of Eq. (El). 
First, the validity of Eq. ( ]A2I) is independent on an approximation for f^c, and it can be 
verified directly with use of the exphcit representation of xks of Eq. (fTTj) . Second, Eq. (1A3P 




(Al) 




(A2) 




(A3) 
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holds if Eqs. (1A1|) and (lA2p hold as can be seen by easily inverting the arguments of Ref. [13 
leading from Eqs. flA2p and flA3l) to Eq. flAl|) . Equation flA3|) can be rewritten as 



j Xik{r,r',uj)V,V^Vo{r')dr' = 

c y Xifc(r,r',u;) Vfer^ dr -no(r)5ij. (A4) 

The next step involves integration by parts requiring the response function to vanish at infin- 
ity and, therefore, it applies to bounded systems only. In this case we can write multiplying 
Eq. (IA4p scalarly from the left by V and using Eq. ([3D 

j x{r,r',u) [u% - V'^Voir')] dr' = V,no(r), 

and after the inversion 

j x"'(r, r', cu)V,rio(r') dr' = cuh, - V,Vo(r). (A5) 

A similar relation holds for xks 

J XkU^, r', Lu)\/M^') dr' = - V^VKsi^). (A6) 



Subtracting Eq. flA5l) from Eq. flA6l) and using the definition of Eq. ([T]), we immediately 
arrive at Eq. (j6]). 
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